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ABSTRACT 

Using large number of Rossi X-ray Time Explorer observations of Sco X-1 we present a detailed investigation 
of the transition layer (TL) and the relativistic precession (RP) models. These models predict the existence of the 
invariant quantities: an inclination angle 8 of the magnetospheric axis with the normal to the disk for the TLM 
and a neutron star (NS) mass M^^ for the RPM. Theoretical predictions of both models are tested and their self- 
consistency is checked. We establish that: (i) The inferred 5— value is 5°.56 ±0°.09. Correlation of the 5— values 
with the horizontal branch oscillation (HBO) frequency is rather weak, (ii) There is a strong correlation between 
an inferred and the HBO frequency in the RPM frameworks, (iii) We infer for different assumptions 
regarding the relations between the HBO frequency V^g^ and the nodal frequency V^^^^. We find that the infeiTed 
Mj^^ = (2.7 ± 0. 1 )Mp, cannot be consistent with any EOS of NS matter We conclude that RPM fails to describe 
the data while TLM seems to be compatible. 

Subject headings: Accretion, accretion disks ■ 

1. INTRODUCTION 



■ relativity — stars: individual (Sco X-1) — stars: neutron 



Kilohertz quasi-periodic oscillations (QPOs) have been dis- 
covered by the Rossi X-ray Timing Explorer (RXTE) in a num- 
ber of low mass X-ray binaries (Strohmayer et al. 1996, van 
der Klis et al. 1996). The presence of two observed peaks with 
frequencies Vj and in the upper part of the power spectrum 
became a natural starting point in modeling the phenomena. At- 
tempts have been made to relate Vj and and the peak sepai'a- 
tion Av = V2 — Vj with the neutron star (NS) spin. In the beat 
frequency model (Miller, Lamb & Psaltis 1998) the kHz peak 
separation Av is considered to be close to the NS spin frequency 
and thus Av is predicted to be constant. In other words Av is an 
"invariant" which does not vary when the kHz QPOs frequen- 
cies change. However observations of kHz QPOs in a number 
of binaries (Sco X-1, 4U 1728-34, 4U 1608-52, 4U 1702-429 
and etc) show that the peak separation decreases systematically 
when kHz frequencies increases (see a review by van der Klis 
2000, hereafter VDK). For Sco X-1 VDK found that the peak 
separation of kHz QPO frequencies changes from 320 Hz to 
220 Hz when the low peak Vj changes from 500 Hz to 850 Hz. 
The lower frequency part of the power spectrum, contains two 
horizontal branch oscillation (HBO) frequencies V^^^ ~ 45 Hz 
and V2HB0 ^ (presumably the second harmonic of V^g^) 

which slowly increase with increase of Vj and Vj (van der Klis 
et al. 2000). Any consistent model faces a challenging task 
of describing the dependences of the peak separation Av and 
HBO frequency VfjgQ on Vj and V2 . 

1.1. QPO model description 

There are two other QPO models in the literature which infer 
the relations between Vj, Vj and V^^^. 

The transition layer model (TLM) was introduced by 
Titarchuk, Lapidus & Muslimov (1998), hereafter TLM98, to 
explain the dynamical adjustment of a Keplerian disk to the 
innermost sub-Keplerian boundary conditions (e.g. at the NS 
surface). TLM98 argued that a shock should occur where the 



gested by Osherovich & Titaixhuk (1999), hereafter OT99, that 
the radial oscillations of the fluid element that bounced from the 
disk shock region (presumably at the adjustment radius) would 
be seen as two independent oscillations in the radial and the 
vertical directions due to the presence of a Coriolis force in 
the magnetospheric rotational frame of reference. Simultane- 
ous measurements of the frequencies of the two kHz QPOs and 
HBO harmonics in a wide frequency range allow one to derive 
the angle 5. Titarchuk & Osherovich (2001, hereafter TOOl), 
claim that for the sources GX 340H-0, Sco X-1, 4U 1702-42 and 
4U 0614H-09 the inferred angle 5, (see Eq 1 in TOOl) stays the 
same over significant range of the observed QPO frequencies. 
The low branch frequency v^, the Kepler frequency Vj^ and the 
hybrid frequency v^^, as they are introduced by OT99 are eigen- 
frequencies of the oscillator However these frequencies are re- 
vealed in the observations as the resonance frequencies V^^^, 
Vj and V2 which are broadened as a result of the (radiative) 
damping in the oscillator (see TLM98, Eq. 15). Furthermore, 
the resonance frequencies V^g^, Vj, Vj are shifted with respect 
to the eigenfrequencies v^, Vj^, V^^. The frequency shift and ran- 
dom errors of the eigenfrequencies depend on the damping rate 
of oscillations A (see details in §2). One should keep in mind 
that the systematic and random errors in the centroid frequency 
determination due to this resonance shift can be a factor of a 
few larger than the statistical error in the determination of the 
centroid frequency. In the present data analysis (as the first ap- 
proximation) we assume that Vj^ = Vj , 
We discuss the results of this analysis in §2. 

The relativistic precession model (RPM) is related to high 
speed particle motion in strong gravitational fields, leading to 
oscillations of the particle orbits. Bardeen, Press & Teukolsky 
(1972), Okazaki, Kato & Fukue (1987), Kato (1990) and later 
Stella et al. (1999) studied precession of the particle orbit under 
influence of a strong gravity due to the GR effects. In order to 
generate these oscillations one should assume that the particle 
orbit is not precisely in the equatorial plane of the compact ob- 
ject. The inclination angle between the particle orbit plane and 



V,, = V2 and v^ = VfjBo- 



Keplerian disk adjusts to the sub-Keplerian flow. It was sug- 
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the equatorial plane can be infinitesimal. For a nonrotating con- 
figuration of the central body (in the Schwarzchild's treatment) 
the angular velocity is calculated using a classical Kepler's for- 
mula. In the general case, when the intrinsic angular momen- 
tum, fl 7^ the azimuthal frequency is expressed in units G = 1, 
c = 1 as follows (Bardeen et al., 1972; Stella et al., 1999) 

V^ = {M/ry/^[27t[l+aiM/ry/']y\ (1) 

Hereafter, we consider a corotating neutron star and a probe 
particle, namely with v,^ > and a > 0. An epicycUc fre- 
quency Vr along with determine the orbital periastron ro- 
tation Vper = Vp — Vr. A preccssion frequency of the node line 
^nod calculated using and a frequency of oscillations in 
the direction which is perpendicular to the orbital plane, Vg, 
namely V^^^ = — Vg|. Formulas for Vr and Vg: 

Vg = V^[l -4fl(M/r3)l/2 + 3fl2/^2jl/2^ 

-8fl(M/r^)i/2„3^2/^2ji/2 



see details in Vihlinin et al. 1994) and a few Lorentzians to 
determine QPO frequencies and we added the following term 
P{v) =Av""[l + (v/Vi^^^^i^)^]^^ , to our model to approximate 
low-frequency noise component. 



5.5 



nur = vJl ~6M/r- 



(2) 
(3) 
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. — Dependence of the inferred 5 angle on the HBO frequency, 
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have been obtained by Okazaki, Kato, & Fukue (1987) and 
Kato (1990). In the framework of the RPM (e.g. Stella et 
al., 1999) a frequency of Keplerian rotation is related to the 
higher kHz peak Vj, and a frequency of the periastron preces- 
sion Vper = Vq - Vr is related to a lower kHz peak Vj ^. The 
HBO frequency V^^^ is related to the nodal precession V^^^^. 
Stella et al. (1999) proposed that v^^^j can be an even har- 
monic of v„^j. Thus equations (l)-(2) allow one to find the 
neutron star mass Mj^^ and a using the observable frequencies 

The goal of this Letter is to test the invariant predictions 
of the TLM and the RPM using extensive RXTE observations 
of QPO phenomena. The Letter reports the results of the de- 
tailed data analysis from Sco X-1 collected by RXTE during 4 
years of observations from 1996 to 1999. In §2, we describe 
the RXTE data which we use to construct the power spectra in 
the frequency range from ~ 0.03 Hz to ~ 2050 Hz. In §2 we 
also give details of the resonance effect on the eigenfrequency 
restoration using the observed QPO frequencies and present 
comparisons of the predictions of the QPO models with the 
RXTE observations. Sunraiary and conclusions are drawn in 
§3. 

2. OBSERVATIONS AND DATA ANALYSIS 

We used data of from the RXTE proportional counter array 
(PC A) instrument (Jahoda et al. 1996) retrieved from the high 
energy astrophysics archive of NASA/GSFC. Sco X-1 was ob- 
served by RXTE during 11 observations (10056, 10057, 10059, 
10061, 20053, 20426, 30035, 30036, 30406, 40020, 40706) 
in 1996-1999. Data were collected from either energy chan- 
nels 0-87 (-2-20 keV) or 0-249 (-2-60 keV) at resolu- 
tions of — 244/is or higher. For observations when "double- 
event" data^ were available, we combined them with single- 
event to increase statistics. We constructed the power den- 
sity spectra (PDS) in the frequency range from 0.03125 Hz to 
2048 Hz to analyze the Sco X-1 variability. The features of 
the low frequency and high frequency variability were inves- 
tigated in the range 0.1 — 256 Hz and 256 — 2048 Hz respec- 
tively. The PDS was fit by a constant (due to dead time effect. 

In the Kerr's geometry ^ Vg (for the Schwarzchild's case with a dimensionless angular momentum a = these frequencies coincide = v^). 

' Two events detected within 6fts at two different anodes of detector are due to charged particles and high count rate from very bright sources like Sco X-1. In this 
case the part of the source flux dominates in double-events data. 



ks. Upper panel: that presented in upper panel is averaged over 5 HBO fre- 
quencies. Approximations by a constant value are shown by solid lines. 



2.1. HBO frequencies and its harmonics 

Among all analyzed spectra where a kHz QPO pair was de- 
tected we selected only those for which an HBO was detected 
as well and at a significance more than 4 G. The second har- 
monic VjfjgQ is present for nearly all of these V^^^. It is worth 
noting that the significance of the HBO peaks decrease with the 
frequency which is a well known effect for quite a few sources 
noted by VDK. Because of this, we cannot find HBO harmonics 
for HBO frequencies higher than 48 Hz (Vj^^^ > 96 Hz) with 
reUable level of significance. Ratio fj^j = '^2HBoI^hbo' 
average, is slightly below than 2, namely ''2/1 = 1 .965 ± 0.004. 
The observational appearance of the HBO harmonics is pre- 
sumably a combination frequency effect along with a reso- 
nance in weakly nonlinear oscillation system (Landau & Lif- 
shitz 1965, hereafter LL). 

2.2. Resonance effect 

The resonance in weakly nonlinear systems occurs at the 
eigenfrequencies of the system Vq when the frequency of the 
driving force 7 is 7 « P'^^/q and p, q are integers (LL). The 
main peak resonance power (for p = 1 and q = 1) is the 
strongest among all the harmonics because the peak power di- 
minishes very quickly with the increase of p and q (LL). The 
maximum of the main resonance amplitude for the linear and 
nonlinear oscillations is not precisely at the eigenfrequency Vq, 
but rather it is shifted to the frequency Vr = Vq -|- 7] . For a Unear 
oscillator Vr depends on the damping rate X (see LL and and 
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TLM98, Eq. 15 for the power of the forced oscillations), i.e. 
Vr = (Vq - 2A^)^/^ and thus the frequency shift 

tj = (Vo2-2A2)V2_Vo. (4) 

From Eq. (3) it follows that ?] « -A^/Vg = -A/g when 
2 = Vq/A > 1. LL also shows that the half- width of the 
resonance curve Av^^y w A for g ;» 1. Taking a differ- 
ential from the left and right hand sides of the formula 
for Vr one can derive a formula for the random error of 
the eigenfrequency determination AVq due to the resonance 
effect: AVq = VrAVr/Vg -h 2AAA/Vq. Now one can esti- 
mate the ratio of the systematic shift and random error to 
^^QPO- They are TT/AVg^^ « Av^^/QAVg^^-S - 10 and 
^^oMvq/.o > 2Vf5yAVf^/(VgpoAVgpo)>2-3 for frequen- 
cies presented by Jonker et al. (2002) in Tables 2 and 4 ^. 
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(^HBO + l^lHBo)- panel: 5 for a particular observation 
with a duration ~3 ks. Upper panel: that presented in upper panel is averaged 
over 5 HBO frequencies. The best fits of the 5— values are shown by solid 
lines. 

We can conclude that Jonker's et al. for the 5— mean would 

be significantly reduced if this effect were taken into account. 

In the nonlinear case the systematic and random shifts are 
also affected by the amplitude of the driving force and the 
relative weights (coefficients) of the nonlinear terms (LL). It 
is not by chance that our observed values of r2^j differ from 
2, since the apparent observed resonance frequency V^^q is 
shifted with respect to the eigenfrequency Vq. This realiza- 
tion is particularly important when one compares the observa- 
tionally identified frequencies V^^gQ, Vj and Vj and those pre- 
dicted by the QPO models. In order to lessen the uncertainty 
in determining the real eigenfrequency value as well as to de- 
crease the systematic error caused by this nonlinear effect, we 
use two values for the low-frequency identification: Vfj^Q and 



'HBO 



calculated using the main resonance frequency v^j^q and the 
second harmonic frequency V2hbo- average PDS and the 
main frequencies v^^^j (and also v^^o^' ^2 ^^^^ used 



for the analysis of the TLM and RPM. 

2.3 . Transition Layer Model (TLM) 
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Fig. 3. — Dependence of the inferred neutron star mass, M^^j on the HBO 
frequency, v^g^. The used data are the same as for Fig. 1 Approximations by 
a constant value are shown by solid lines 

For each set of v^^^, Vj and we obtained the 5— value us- 
ing formula (1) in TOOl, where VfjgQ, Vj and Vj are identified 
with v^, V]^ and v^ respectively. In Figure 1 the dependence of 
the angle 5 between a normal to the disk and a magnetospheric 
axis on V^^^ is shown. We establish that the strongest depen- 
dence of 5 on the QPO peak frequencies involves Vfj^g- Thus 
we investigated this dependence only. The model invariants 
should be kept the same, independent of the QPO frequency 
branch where it is verified. In Table 1 the approximation of 5 
by constant are present for two cases: v^ = VfjgQ (see Fig. 1 the 

lower panel) and v^ = v^^q = \ ("^hbo + 5 ^2hbo) ^^^^^^ 
case the second harmonic is taken into account for 5 determi- 
nation). For each case we provide a reduced Xyg^- The 5-^rms 
and the standard deviations Gg are also present in Table 1 . We 
repeated the same procedure for 5 vs v^ when v^ = v^^q. In 
Figure 1 these results are present. In order to reduce the pos- 
sible stochastic scattering in the data, which always dominates 
relative to the statistical scattering, we implement an averaging 
of 5 for a given HBO frequency (see the upper panel of Fig. 1). 
We repeated the same procedure for 5 vs v^ when v^ = v^^o- 
In Figure 2 these results are present. Whereas some correlation 
of 5 vs VfjgQ are clearly seen in Figure 1, points (5, v^j^q) are 
almost scattered along the horizontal line in Figure 2 (see the 
upper panel). 

2.4. Relativistic Precession model (RPM) 
Now we investigate the RPM invariant, the NS mass 



HBO 



\od or as a function of V^^^ 



nod' 



The deduced dependences of M^^ on the HBO frequencies 
are stronger than for the kHz frequencies (similar to the case 
with 5 for the TLM). In Table 2 and in Figs. 3-4 we present 



* Jonker et al (2002) presented the analysis of the 5 - invariant of the TLM. They found that 5 = 6" .3 ± 0°. 1 and 5 = 6°. 1 ± 0°. .2 for GX 5- 1 and GX 17+2 respec- 
tively. But they pointed out the very large Xred ~ 95.2/9 and ~ 377./9 for GX 5-1 and GX 17+2 respectively. Their calculations include the statistical error 
only but do not account for the systematic resonance shift T] and for the random resonance error AVg. 
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a mean value of M^^, Xred ^^'^ "m- ^^'^^ Xred values 

are a factor of 3-6 bigger than that in Table 1 for 5. But the 



relative deviation of M, 



is also 2% as Og for 5. The lin- 



ear correlation of the Mj^^ values on HBO frequencies v^gg is 
clearly seen in Fig. 3 and gets even stronger for V^^^ (Fig. 4). 
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Fig. 4. — Dependence of the inferred neutron star mass on the HBO fre- 
quency, v^g(j. The used data are the same as for Fig. 2. Approximations by a 
constant and linear trend are shown by dashed and soUd lines respectively 

increase 



NS 



HBO 



of V 
mean M, 
in the 



drops from 2.8 
from 39.5 Hz 



Mg to 2.6 M, 
to 46.5 Hz 



© 



7 



with 
In RPM 



the 



are scattered along the oblique line whereas 
TLM 5 are scattered along the horizontal Une. 

Table 1 



Parameter 


V,, = v„«„ 


V,, = v„«^ 


S 


5.612±0.007 


5.545±0.006 


^red 


5.53 


2.94 


s 


5.600±0.114 


5.555±0.086 




0.567 


0.359 


Proh 


210-'' 


310-2 


N, bins 


39 


39 




0.738 


0.595 


Prob 


3.7-10-2 


1.2-10-1 



N, bins 



See caption for Figs. 1-2: v^gQ - 
eter, 5= l^S,.. 



-l-iv„ 



5 ), 5 is the best-fit param- 



2.5. Statistical evaluation of the suggested invariants for the 
TIM and the RPM 

In order to check the consistency of the TLM and the RPM 
we use the Spearman's non-parametric test for evaluation of 
the correlation of 5 and M^^g with HBO frequencies. If the 
correlation coefficients = ±1, then there is a perfect linear 



correlation between two values. In the general case |rj| < 1. 
For each case we also calculate a probability Prob that the 
derived correlation is a result of the statistical fluctuations ^. 

Table 2 



Parameter 








2.680±0.002 


2.675±0.002 


^red 


14.3 


16.3 




2.677±0.050 


2.674±0.053 


a, cm X 10^ 


1.367± 0.002 


1.358±0.002 




8.2 


13.7 


N, bins 


39 


39 




-0.715 


-0.789 


Prob 


3- 10-7 


2-10-9 




-0.929 


-1.000 


Prob 


8.6- lO-"* 


0.0 


N, bins 


8 


8 



See caption for Figs. 3-4: a = J/{c ■ M^j), 
is best-fit parameter, M^,^- = jj Y,Mns , - 



^y^HBO'^ 1^2HB0 



From Table 2 it is seen that the RPM is inconsistent with data 
(r^ = — 1, there is a perfect linear correlation between M^^ and 

) and the proposed parameter M^^ is not invariant with 



'HBO 

respect to the observed QPO frequencies . It is worth noting that 
inclusion of the second HBO component into the determination 
of the HBO peak leads to a weakening of the correlation for 
5 vs HBO frequencies while the opposite effect is seen for the 
correlation between vs HBO frequencies. In Figs. 2, 4 it is 
shown that the averaging leads to a decrease in the scattering of 
the inferred parameters along the best-fit line. For the TLM (see 
Fig. 2) 5 is practicafly a constant, 8 = 5°.55 with « 0.5% 
while for the RPM (Fig. 4) strongly correlates with V^j^q 
and is inconsistent with a constant. It should be pointed out that 
in the Table I are affected by the aforementioned under- 
estimation of the systematic and random resonance shift (see 
Introduction and §2.2). Thus, the j^^^^ is artificially high be- 
cause it takes into account only statistical errors, and ignores 
this effect. 



3. CONCLUSIONS 

Using the thorough data analysis of RXTE data for Sco X-1 
we have presented a detailed investigation of the transition layer 
(TL) and the relativistic precession (RP) models. These models 
predict the existence of the invariants: 5 for the TLM and M^^ 
for the RPM. We establish that: (1) The inferred 5— value is 
5°.56 ± 0°.09 is consistent with a constant. Correlation of the 
5— values with the HBO frequency is rather weak. (2) There is 
a strong correlation between an inferred neutron star mass and 
HBO oscillation frequency in the RPM frameworks, (3) The in- 
ferred NS massM^^ = (2.7±0.1)Mq cannot be consistent with 
any EOS of NS matter. 

We are grateful the referee for the useful suggestions and 
the valuable corrections. L.T. acknowledges the fruitful dis- 
cussions with Chris Shrader, CharUe Bradshaw, Vladimir 
Krasnopolski and Paul Ray. 



7 As it is expected from the probability theory, the statistical mean Mj^^ and the statistical deviation converge to the real values of the mean and the standard deviation 
with increase of a number of points for a given HBO frequency. 

* Particularly, this test is useful for comparison whether one of the correlations is stronger than another. 
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